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Fully twenty years ago Grady et al. (1987) and 
Sonneborn et al. (1987) reported that the UV reso- 
nance lines of the sharp-lined B2 IV-V star £ Cassiopieae 
(HD 3360;V=3.67) are variable. Interest in this star has 
been rekindled recently by the report by Neiner et al. (2003; 
hereafter N03) that this star has a dipolar magnetic field of 
B 2 w 335 G varying by ±28 G with a period of 5.37 days. 
The resonance lines of N V, C IV, and Si IV show variations 
with the same period. N03 also report a nitrogen abun- 
dance enhancement and a marginal helium overabundance 
([He] = 0.11 ±0.06 dex). These characteristics make CCas 
a likely weak member of the magnetic He-strong Bp stars. 
We note for completeness that small-amplitude nonradial 
pulsations are excited on this semi-evolved star's surface 
(e.g., Sadsaoud et al. 1994). 

It is generally accepted that magnetic Bp stars con- 
tain fossil dipolar magnetic fields that are randomly ori- 
ented with respect to the rotational axes. Many authors 
(e.g., Shore 1987, Shore & Brown 1990, Babel & Montmerle 
1997a, 1997b) have attributed the time-modulated UV res- 
onance line variations to the channeling of radiative winds 
from their magnetic poles to their equators. Because ionized 
particles cannot cross magnetic field lines, the global stel- 
lar fields play a key role in the wind outflows. Although the 
particles emerging from the magnetic poles can freely ac- 
celerate, wind flows are prevented from crossing transverse 
stellar field lines in the equatorial zones and are suppressed. 
Moreover, the field lines emanating from the stellar mag- 
netic poles channel the escaping wind toward the magnetic 
plane. Here they collide with wind streams coming from 
the opposite magnetic hemisphere. This produces a shock 
that overionizes particles which then reionize and deexcite 
and emit resonance line radiation. The resonance line ra- 
diation is then scattered through the confined structure, 
allowing the observer to see an array of absorption and 
emission components as the corotating disk/shock complex 
corotates around the star. Many of the cooled particles set- 
tle toward the plane to form a confined "disk." Interior to 
the magnetospheric structure, a sector of this disk transits 
in front of the star and occults part of the star's projected 
area and causes excess absorption in low-excitation lines of 
iron group ions. In Bp stars with extensive disks, the ab- 
sorption from this "iron curtain" forest is made visible by 
their composite effect in high resolution spectra obtained 
by the International Ultraviolet Explorer (IUE) (Smith & 
Groote 2001; hereafter SG01) and in high-level Balmcr lines 
(Groote & Hunger 1976). 

In the hydrodynamic formulation of this process by 
Babel and coworkers, a shock develops where the disk par- 
ticle pressure equals the wind's energy density. Important 
refinements have been added to the Babel picture by ud- 
Doula and collaborators by considering the relative influ- 
ences of other forces on particles along the flow. We sum- 
marize several effects they can have on the disk structure: 

- Ud-Doula & Owocki (2002; hereafter uO02) found that 
the confinement of the wind by the star's magnetic field 
is given by a ratio "77*" of surface magnetic to wind en- 
ergy densities. These authors also found that the weight 
of the wind particles produces a feedback on the field 
lines in which the stagnated wind drags the lines of 
force inward along the equatorial plane. In some circum- 
stances, the particles may collapse the field lines back 
into the star's equatorial regions. More typically, con- 



densations form chaotically as snake-like blobs. These 
structures fall back toward the star, out of the plane 
and along the lines of force. UO02 suggested that this 
sequence of events is responsible for the UV emission 
lines in Bp star spectra. Gagne et al. (2005; G05) devel- 
oped this idea to explain emission from the hot Bp ana- 
log 8 1 OriC. A byproduct of this picture is that matter 
does not fall straight inward along the magnetic plane 
to the star's surface. The effect is to produce a discrete 
inner edge to the disk. Interior to this point there is 
effectively no matter. 

— Refinements in the picture have shown the importance 
of angular momentum conservation and the tendency 
of particles to accumulate in magnetospheric regions of 
minimum potential energy (Preuss 2004, Townsend & 
Owocki 2005). These forces can cause the disk to warp 
for a range of magnetic inclinations and also permit the 
particles to accumulate at the intersections of the mag- 
netic and rotational equators and in other areas outside 
the planes. The effect of these accumulations is that 
disks of rapidly rotating Bp stars may show departures 
from coplanarity and axisymmetry. 

— Ud-Doula et al. (2006; u06) found that lines of force near 
the Alfven point of rotating stars are stretched outward 
by the centrifugal force. At irregular intervals, the field 
lines severe, and disk material "breaks out" through the 
Alfven point at the equator. The heating of material by 
this process may be the source of X-ray flares, such as 
observed in the prototypical magnetic Bp star a OriE 
(e.g., Groote & Schmitt 2004). 

— Including post-shock cooling energetics into quasi- 
equilibrium magnetic wind models leads to the predic- 
tion of a range of temperatures from 10 8 K to 10 6 K or 
less (Tonneson et al. 2002, G05). This shock region is 
confined to a broad zone just inside the Alfven radius. 
Interior to this zone, wind particles are not accelerated 
to high enough velocities to produce "superions" and 
resonance line radiation. 

For relatively slowly rotating, weakly magnetic stars like 
C Cas some of these processes are not expected to be impor- 
tant. An important modification of the Babel model is that 
the wind detritus returns to the star along lines of force 
where it collides with high latitude wind streams. 

1. Methodology 

1.1. Reduction of ultraviolet resonance lines 

The ultraviolet spectra for this program included the wave- 
length regions surrounding the N V, C IV, Si IV, and Al III 
resonance lines and recorded on IUE echellograms with the 
Short Wavelength Prime (SWP) camera. These data were 
obtained from the MAST archivesQ We found 106 spec- 
tra, all obtained between 1979 and 1995 and all observed 
through the IUE's large aperture, were suitable for analy- 
sis. We then computed "equivalent width (EW) indices" 
by ratioing the fluxes within the line edges with fluxes 
of other wavelengths in the order (see SG01). Spectra for 
phase ranges given in tJ2.1l were selected, "conditioned," and 
coadded by carrying out the following operations: (1) cross- 
correlating them to remove apparent wavelength shifts, (2) 
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forcing them to a common continuum slope across the or- 
der, and (3) weighting the individual spectra according to 
the mean fluctuation from pixels of neighboring resolution 
elements. This weighting was accomplished by evaluating 
the internal spectrum r.m.s. noise (assumed to be the me- 
dian value of the absolute value of the fluctuations of pix- 
els in each neighboring resolution element, i.e., every 3rd 
pixel). The principal work in this paper concerns the anal- 
ysis of pairs of coadded resonance line spectra obtained at 
phases near the passage of the South magnetic pole across 
the star's central meridian (following N03, taken to be <j> 
= 0.0), and approximately 0.5 cycles later, i.e., when these 
lines exhibited minimum and maximum absorption, respec- 
tively. We quantified the variations of these line profiles in 
our first figure by measuring a dimensionless "equivalent 
width index" rather than a true equivalent width (EW) in 
Angstrom units. Our index is defined as the ratio of all 
fluxes within a designated window across the profile to all 
fluxes outside them in the echelle order containing them. 
Our indices are proportional to a true equivalent width but 
do not rely upon inaccurate placements of the local con- 
tinua. For the case of £ Cas, the windows were taken to be 
the instrumcntally broadened photospheric line profile, or 
±0.8 A for most resonance lines (±0.5 A for N V). 

1.2. Optical Data 

In an attempt to search for variations in high-level Balmcr 
lines from the variable occultation of the star by the 
disk, we obtained 31 spectra with the 1.85 m telescope 
of the Dominion Astrophysical Observatory during 2005, 
December 16, 2006 February 14-16 and March 8. We used 
the f/5 Cassegrain spectrograph and fed the light to a 
CCD detector by means of an image sheer. The result- 
ing spectra have a dispersion of 10 A mm -1 , a resolution 
of about 16 400, and cover the wavelength region AA3610- 
3880. These data gave a null result, as reported below. 

1.3. Spectral Synthesis Computations 

To perform the analysis of the IUE spectra, we utilized 
a suite of programs written by I. Hubeny and collabora- 
tors and non-LTE solar abundance BSTAR models (Lanz 
& Hubeny 2006). The first of these, synspec, is a photo- 
spheric line synthesis program (Hubeny, Lanz, & Jeffery 
1994) used for the reference photospheric profiles. This pro- 
gram permits one to modify assumed chemical abundances 
in the atmosphere. Indeed, the strengths of the Felll lines 
in the IUE spectrum suggested an iron abundance of ^ 
the solar value, assuming a photospheric microturbulcnce 
£t = 2 kms -1 . We also used synspec to convolve the photo- 
spheric lines with functions approximating the instrumental 
and rotational broadenings for the instrument and the star. 
The spectroscopic resolution of IUE high-dispersion spec- 
tra is about 13 000 and the the star's rotational broadening 
(vsini = 17 kms -1 ) arc comparable . We computed the 
synthetic spectra in steps of O.OlA. 

Our goal in this analysis is to assess the effects of the 
magnctosphcric disk on the UV resonance line spectrum as 
the disk moves across the surface of ( Cas, alternately oc- 
culting greater or lesser areas of the star in its circuit. To 
simulate the absorption and emission signatures of the disk 
on the spectrum, we used the radiative transfer program 



circus (Hubeny & Heap 1996, Hubeny & Lanz 1996). This 
program computes strengths of absorption or emission com- 
ponents of lines in a circumstellar medium from user-input 
quantities such as disk temperature, column density, areal 
coverage factor (the portion of the disk that passes in front 
of the star), and microturbulent velocity. Circus combines 
this contribution with the sj/nspec-computed spectrum of 
the photosphere. In its solution of the radiative transfer, 
circus calculates line emissions and absorptions separately 
along a line of sight. The emission is computed according 
to the assumption of LTE for an input disk temperature 
Tdisk- Because the "absorption" of resonance lines is actu- 
ally due to scattering, it can be handled by turning off a 
flag in the program that computes the emissivity coefficient 
in the radiative transport equation. The result is that the 
background flux of the star is affected only by the absorp- 
tion coefficient. The ability of handling the line transfer in 
this way permits one to treat absorptions and emissions si- 
multaneously or separately in different circumstellar struc- 
tures. Circus can accommodate three separate "clouds." In 
the present analysis, we found it necessary to include two 
cloud structures, one for the excess absorptions of the reso- 
nance lines and a second for the emission components. The 
program allows us to input different radial velocities and 
temperatures that might be useful in testing whether dif- 
ferent diagnostics are formed in separate regions, allowing 
an assessment of whether they might arise in different vol- 
umes. As in the work described by SG01 and Smith et al. 
(2006; hereafter S06), we have assumed the clouds to be 
homogeneous in all properties. However, the MHD studies 
described above suggest this may well be an oversimplifica- 
tion. 

Among other input parameters are the disk's volume 
density and microturbulent velocity £. Our results were 
insensitive to the values chosen for either of them. The 
volume density enters only indirectly (and insensitively) 
through the gas ionization equilibrium. Our trial input den- 
sity, ~ 10 12 cm~ 3 , is the order of magnitude SG01 deter- 
mined for hot Bp star disks and is also the peak value deter- 
mined by the MHD simulations. Except for the absorptions 
of the primary resonance lines (C IV, Si IV) the overwhelm- 
ing majority of the UV metal lines in our computations are 
optically thin. Since the spectral signatures then scale as 
the volume of absorbers, a degeneracy arises between two 
geometrical factors of the disk: its area and column length 
along the observer's line of sight. 

The final ingredient of our models is the treatment 
of line of sight velocities. Because the magnetosphere co- 
rotates over the star's surface, no differential velocities 
exist along the line of sight to the star during occulta- 
tion. Moreover, the Doppler shift of the disk absorbers is 
matched by the velocity of the star area behind them. SG01 
found that turbulent velocities in the range of 0-20 kms -1 
can be offset by changes in the assumed metallicity in the 
disk because the lines are optically thin. These authors also 
found a broadening of the Al III and NV resonance lines 
of about 25 kms -1 and 50 kms -1 , respectively, at phases 
near maximum occultation. The circus models were com- 
puted with the same underabundant iron abundance as de- 
termined for the photosphere. We used solar abundances of 
C, Si, and Al, but an overabundance [N] = 0.4 dex (N03), 
all with a £ = 20 kms" 1 . 
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1.4. Relevant rotational and magnetic parameters 

In our analysis we have assumed the physical parameters of 
£ Cas determined by N03. These include the redetermina- 
tion of the star's effective temperature, near 21 000 K and a 
logarithmic surface gravity of about 3.8. Therefore, we have 
chosen the (21,000, 4) model in the BSTAR grid as input to 
our synspec computations. At the surface of the star, small 
amplitude nonradial pulsations are present (N03), but we 
ignore them in this discussion. Given a period of 5.37 days 
and the star's v c sini = 17 kms , the inferred equatorial 
rotational velocity, v c = 55 kms -1 , is also well determined: 
C Cas is a moderately rotating B star. N03 give a rotational 
obliquity angle of i = 18 ±4°. From i and attributes of 
the magnetic phase, curve N03 derive a magnetic obliquity 
angle, (3 = 79 ±3°. 

The rotational phase zero cj) is taken as the stellar lon- 
gitude containing the negative poles in the visible hemi- 
sphere. Then the viewing angle a can be written as: 

cos ot = sin (3 cos <p sini + cos [3 cos i, (1) 

one sees that for small i and large (3 that the viewing angle 
changes little through the rotation cycle. Consequently, we 
can expect only small changes in the absorption coverage 
factors of the disk in front of the star. 

The relation for the wind magnetic confinement param- 
eter 7]*, is given by uO05 as the ratio of magnetic to kinetic 
energies, namely as: 

r h = B 2 R 2 /(M Voo ) . (2) 

Here the symbols refer the the polar magnetic flux B z , stel- 
lar radius R*, and mass loss rate M, and the wind velocity 
at infinity Voo. One is somewhat handicapped in evaluat- 
ing 77* because no modeling has been performed of C Cas's 
wind attributes. Therefore, we will take mean values for a 
B2V star: M = 3xl0~ 10 Moyr -1 (Snow 1981). We es- 
timate that = 800 kms -1 from the line profiles and 
i?* = 5.9R Q from N03. Then, 77* -10 3 . This is a large 
value but still about 10~ 4 times the value found for a OriE 
and 36Lyn (S06) However, the mass loss rate for £Cas 
may be anomalously high. For example, Slettebak (1994) 
found that the combined equivalent widths for the the C IV 
and Si IV doublets for this star are 2.40 A and 3.59 A, re- 
spectively, whereas the corresponding values for an average 
B2 V star are 1.35 A and 2.85 A. If this is the case, then 77* 
should be taken as an upper limit. Finally, we can estimate 
the location of the Alfven radius from the approximation 
Raij ~ rj*. Thus, RAif,* ~ 5-R*. 

2. Results 

2.1. The UV resonance line variations with phase 

N03 have determined a period of 3.57045±0.00008 days 
from the variations of the CIV, Si IV, and NV resonance 
lines. To check this value, we applied the PDM period- 
finding algorithm (Stellingwerf 1978) to EW indices ex- 
tracted from the resonance lines as defined above. These ex- 
tractions were done in the inner ±0.5 A of the N V lines and 
inner ±0.9 A of the CIV and Si IV lines. We ordered these 
values by time and and found a value of 3.57038 ±0.0001 
days, thereby verifying the N03 value. We used this pe- 
riod and the N03 zero point (</> = at HJD 2446879.889) 
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Fig. 1. Variation of line strength indices for the UV res- 
onance doublets of N V, Si IV, and C IV with phase in the 
spectrum of £ Cas. Scaling factors are applied to show the 
same apparent amplitude; vertical shifts are made for clar- 
ity. Approximate measurement errors are indicated. The 
dashed, tapered line shows the assumed limiting absorp- 
tion relation (see text). 

to compute rotational phases. Figure[T] exhibits the absorp- 
tions (EW index) folded around this period, but we have 
rescaled them so that each exhibits the same fractional vari- 
ation (this rendition artificially increases the apparent noise 
in the NV curve). Note the Fig. [T] minimum centered at 
phase 0.0. Also, as previously pointed out by N03, one may 
discern a hint of a secondary minimum on the broader max- 
imum at <j)Ki 0.50, but we are not convinced that it is statis- 
tically significant since this feature does not appear in the 
Si IV or N V curves. Finally, we point out that the shapes of 
the absorption-phase curves are indistinguishable from one 
another. This is surprising because, as noted below, the de- 
creased absorptions they represent occur for both optically 
thick and thin lines. 

To reconnoiter the modulation of the line strengths 
with phase, we first grouped over 20 spectra each into 
EW-minimum and EW-maximum phase bins correspond- 
ing to 0=0.89-0.11 and 0.30-0.70, respectively, and com- 
puted mean spectra from these groups. The results are ex- 
hibited as thick solid lines in Figures^] and for the wave- 
length regions surrounding the CIV and NV doublets, re- 
spectively. Minimum and maximum spectra were also com- 
puted for the Si IV lines. The line profiles of the spectra 
are different in that those for the former phases exhibit 
smaller absorptions and generally steeper red wings. Part 
of the asymmetry in the wings is due to the influence of 
the winds. As for other Bp stars and 6 |1 OriC (Smith & 
Fullerton 2005), this wind is focused by the magnetic field 
and streams outward preferentially at low magnetic lat- 
itudes. Since the wind escapes from regions outside this 
disk, we do not model it in this work. The red shoulders 
of the profiles at these times are steep with respect to our 
computed photospheric profiles, and we attribute this to an 
incipient redshifted emission. We attribute the strengthen- 
ing of the spectra to increased blocking of the stellar flux 
by the intervening disk. Our modeling of absorptions and 
emissions proceeded with these assumptions, and we defer 
the discussion of results to the next section. 
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2.2. Quantitative analysis of the C IV and N V resonance lines 

2.2.1. The maximum and minimum absorption spectra 

We have described a picture in which wind particles escap- 
ing the magnetic poles of a Bp star shock when they im- 
pact particles near the magnetic equator that have formerly 
settled to form a cooled "disk." The phase-dependent ab- 
sorptions of key ultraviolet lines validates this basic picture 
for C Cas. The morphology of the resonance lines shown is 
replicated within the errors by less excited resonance lines 
of Al III and C II, thereby providing evidence that the struc- 
ture is a confined cool disk. However, we were not able to 
estimate the disk temperature from ratios of iron line ab- 
sorptions because the available Fell lines are too weak to 
show variations. Even so, if we assume the relation for the 
cool disk ionization temperature Tdisk ~ 0.6T o ff, as found 
in the Bp stars studied by SG01 and S06, we obtain a value 
Tdisk — 12 000 K. Assuming this temperature and, in the 
absence of information about the extent of the coverage of 
the star by the cool disk component during the occultation 
phase, wc find from circus modeling of the small variations 
of the A1III A1855-62 and the CII AA1335-6 multiplets 
that the hydrogen column densities of the occulting disk 
are lxlO 18 cm" 2 , assuming £ = 20 kms" 1 , and 2xl0 18 
cm" 2 , assuming £ = 10 kms -1 . We take these as conve- 
nient lower and upper estimates. 

We note that after reducing our DAO spectra of the 
high-level (H8-H13) Balmer lines and plotting the equiva- 
lent widths with phase, we were unable to find statistically 
significant variations in them. From our the measurement 
errors, ±0.6% per line, this undetectability translates to a 
column density of ~ 4 x 10 20 cm -2 , according to our circus 
simulations, and assuming a Tdisk = 12 000 K. Since this 
limit is a few hundred times higher than the limits imposed 
by the All III and C II variations just above, it does not add 
information to the disk analysis. 

The resonance line variations in the two figures bring 
out two interesting characteristics. First, notice that the 
r.m.s. spectra in the C IV line figure show that the profile 
activity is confined to roughly ±100 km s" 1 from line center. 
However, the wind-depressed blue wings undergo no activ- 
ity through the cycle. Similar to the resonance lines of other 
Bp stars, the r.m.s. spectral amplitudes peak to the red of 
the line centers. Inspection shows that this is caused pri- 
marily by steepenings of the red shoulder of the lines in ob- 
servations obtained in the EW- minimum phases. According 
to the line fittings discussed below, these shoulders cannot 
be fit by photospheric spectra or disk absorption spectra. 
Rather, we believe these are incipient emissions correspond- 
ing to the more pronounced emissions observed in many Bp 
star spectra. Another characteristic brought out in these 
plots is that the doublet ratios of these activities are 1.2 
and 1.0±0.2 for the CIV and Si IV (not shown) lines, re- 
spectively. The ratio for the N V lines is 2.5±0.4. From these 
ratios we can see that the circumstellar media in which the 
"disk" component of CIV and Si IV are formed are op- 
tically thick while the NV lines optically thin, much like 
their photospheric components. 

2.2.2. Fitting of the resonance doublets 

In our line profile fitting we found that an important diag- 
nostic is the line optical depths. Because the optical depths 
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Fig. 2. Maximum (upper) and minimum absorption lines 
of the C IV doublet of £ Cas fit with synthetic spec- 
tra (dashed lines). Fitting parameters are indicated (see 
Table 1). The shaded area denotes the emission component 
determined in the models; the dotted line in the lower plot 
repeats the maximum absorption profile in this same wave- 
length range for comparison. The blue wing absorption of 
the C IV lines is due to the (unmodeled) wind contribution. 
The r.m.s. spectra is also shown. 

of our models are typically within an order of magnitude 
of unity, we were able to use them to estimate the optical 
depths of the disk columns and ultimately to set limits on 
the range of the disk temperature Tdisk- These values are 
represented in the table discussed below. 

Our initial attempts to model the minimum and max- 
imum EW spectra in Figs. [2] and [3] used input circus pa- 
rameters for a single homogeneous absorbing and emitting 
circumstellar "cloud." The models assume a disk temper- 
ature Tdisk, which for simplicity we have assumed to be 
the same value. This group of models failed to give sat- 
isfactory fits. Although there were several problems with 
the fits, the most important was the inability of the mod- 
els to produce enough absorption to fit the EW variations. 
This problem was accentuated by the constraint that the 
change of coverage of the star must be incomplete because 
of the low inclination of the star. Investigating this failure 
in more detail, we found that the poor fits extended over 
a significant wavelength range. For this reason, we were 
forced to consider circus models fit with two homogeneous 
clouds with different properties, especially radial velocity. 
These models were successful, but only when we increased 
the microturbulence £ to 75-90 kms -1 . We note by con- 
trast that high cloud macroturbulences merely smear the 
computed excess absorptions and emissions. They do not 
strengthen them, as is required by the observations. 

We were not led immediately to our best fits with the 
two-cloud description. Moreover, those we have selected 
for representation are those obtained only for one possi- 
ble disk temperature among the range of Tdisk values. One 
of the most important fitting parameters, the temperature 
implied by fits to the emission components, T em i s , is not 
uniquely determined. For the C IV line fitting, the upper 
limit on T em i s , 26 000 K, was set by the decreased optical 
depths in the line that led to different and contradictory 
absorbing cloud areas or columns for the two doublet com- 



5 



T 




;V rad (em.)= 75 km/s N v N v 



1236 ~" 1238 " 1240 1242 '" 1244 

Wavelength (A) 

Fig. 3. Maximum (upper) and minimum absorption in the 
N V doublet lines of £ Cas fit with synthetic spectra (dashed 
lines). Fitting parameters are indicated (see also Table 1). 
Vestigial emission (shaded area) in A1242 is suggested in 
our model. 

ponents. The lower limit of 16 000 K was set by the appear- 
ance of nearby Fell lines in the computed models, which 
are not present in the observations. For the N V lines a low 
temperature limit of 26 000 K was again set by the absence 
of a predicted iron line (Felll 1243.4 A). The column den- 
sity required for fitting the NV lines rapidly decreases to 
the values found for the CIV fittings, and we consider this 
the most likely result. 

Table 1 shows the results of best two-cloud models for 
the maximum and minimum EW profiles of the N V, CIV 
and Si IV doublets, that is for observations obtained near 
mid-transit (maximum occultation) phases and those for 
half a cycle later, respectively. An example of a model for 
the CIV and NV doublets is shown in Figs. [5] and [3] as 
a dashed line. Annotations give pertinent model fitting pa- 
rameters. Shadings denote the modeled incipient emissions. 
In the table we give both "high" and "low" Tdisk values, 
as discussed above. In each case we give modeling param- 
eters for fits to the maximum EW and minimum EW line 
in the upper and lower lines of each double row entry. In 
the table we give several key fitting parameters, including 
the Tdisk, column density, areal coverage, radial velocity, 
and central optical depth of the stronger doublet member. 
Errors were estimated by repeating our models with sin- 
gle parameters changed by varying degrees and comparing 
the resulting differences in the equivalent widths with 2<r 
changes in the line profiles within the same wavelength win- 
dow from the r.m.s. errors in the IUE spectra. The latter 
was accomplishing using a Monte Carlo approach to com- 
pute the fluctuations as described in Smith (2006). These 
errors led to the following estimated 2tr errors in our table: 
±15 kms" 1 in microturbulence £, ±20% in fractional areal 
coverage for CIV and Si IV, and ±50% in the column den- 
sities. The errors for the latter two quantities are ±100% 
for the N V doublet. Inspection of this table brings out the 
following points: 

— Unlike the cases of equator-on Bp stars with large 
angles, occultation of some part of the stellar disk occurs 
to one extent or another at all rotational phases. This 



is expected for star-disk systems with a small i. This 
circumstance also hints that the disk is either extended 
in height (perpendicular to the plane) or extends inward 
to reach the star in the magnetic plane. 

— The column density and hence volume ratio of shocked 
gas in which the resonance line emission and scattering 
occurs is several times that of the column density in- 
ferred from the less excited CII and Al III ions. This re- 
sult contrasts with the SG01 finding for other Bp stars, 
where the ratios were comparable. This indicates that 
the ratios of the volumes can be quite different among 
Bp stars and opens the door to the possibility that they 
are not cospatial. 

— An important parameter in our fittings was the radial 
velocity of the emission component - typically ±50- 
125±20 kms^ 1 . Its errors were estimated by a compar- 
ison of values found for the doublet members. 

— The large £ values for Bp disk absorptions and emis- 
sions is typically high (SG01, S06), and we cannot ver- 
ify that it is a true microturbulence. Irrespective of data 
quality limitations, there is an intrinsic interpretational 
ambiguity between geometries that imply a large £ in 
a single structure and multi-columns each with smaller 
£ values and separated in radial velocity. Such a con- 
trivance would both broaden and strengthen the excess 
absorptions and emissions. 

3. Interpretation 

In S06 we described a computer program written to de- 
scribe the varying excess aborptions of an optically thick 
resonance line as a co-rotating disk shaped like a cylinder 
transits across the star. The line is assumed to be optically 
thick, the absorption of flux from the obscured region of 
the star is total. Then, one may express the occulted area 
of the star in terms of the excess absorption measured. 

One can manipulate input parameters representing the 
dimensions of the disk to fit observed phase absorption 
curves. For the case of ( Cas, the curves in Fig.Q] show 
that the absorption curves are essentially flat out to 
±4 cycles around maximum occultation phase, <p ~0.50. 
Unfortunately, flat absorption phase curves give informa- 
tion only on the limits of ratios of dimensions. □ This prob- 
lem is aggravated by the low S/N of IUE data, which makes 
it difficult to map the progress of the disk viewed at a low 
i. However, we may still investigate the general character 
of the disk cylinder shape if we estimate a minimum taper 
in the curves (5% of the maximum excess absorption out to 
4>= 0.25 and 0.75) that might be hidden in the noise of the 
absorption curves. This taper is shown in Fig.[T]as a dashed 
line. 

Applying this criterion and the i and (3 values for £ Cas, 
we adopt an outer disk radius r outer = Ra; / , and we can 
ask, what is the interdependence between the height (mea- 
sured parallel to the rotational axis) and the thickness of 
the disk in the radial direction, r outer - ri nner ? According to 

2 Although we have discounted the statistical significance of 
a secondary dip in the phase curves, if they are real it could 
be interpreted as a dearth of overionized particles like C in 
zones over the rotational poles. This would agree with uO02's 
expectations that the particles exert a feedback on the magnetic 
field and drag them away from zones of low angular momentum. 
Evidence for this avoidance would be visible only for rotating 
magnetic stars observed at low inclinations. 
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Table 1. 



Derived columns, areal coverages, and radial velocities (kms 1 ) inferred from UV resonance lines 



IN V lines 




= 75 km/s): 


















Col. density 




Fractional 


Area 


Vel. 


Vel. 


7~line 








[ADS ) 


[rjimsj 


[ADSJ 


[rjnnsj 


[ADS J 


[Jimis J 


[ADSJ 


[HiimsJ 


34 000 




1.5 x 10 19 




1.0 


0.0 







0.3 




34 000 




1.5 x 10 19 


5 x 10 


0.4 


0.4: 





+75 


0.3 


0.1 


28 000 




3 x 10 21 




1.0 


0.0 







0.6 




28 000 




3 x 10 21 


3 x 10 21 


0.4 


0.4: 





+50 


0.6 


0.6 


Olv lines 




= 90 km/s): 
















Tdisk 




Col. density 




Fractional 


A 

Area 


Vel. 


Vel. 


Tline 


7~Hne 






[ADS ) 


[hjiiiis J 


[ADS ) 


[rjiTiisJ 


[ADS J 


[Jimis J 


(Abs ) 


[Hiims ) 


24 000 




3 x 10 21 




1.0 


0.2 







15 




24 000 




3 x 10 21 


21 

3 x 10 


0.4 


0.8 





+ 125 


5 


5 


on r\r\r\ 

20 (JUL) 




6.0 X 1U 




1.0 


0.0 







40 




ZVJ uuu 




1. / X 1U 


l.i X 1U 


n i 


i n 


n 
u 






n 9 


Si IV A1403 




(C = 75): 
















Tdisk 




Col. density 




Fractional 


Area 


Vel. 


Vel. 


Tline 








(Abs) ' 


(Emis) 


(Abs) 


(Emis) 


(Abs) 


(Emis) 


(Abs) 


(Emis) 


22 000 




5 x 10 18 




1.0 


0.0 







2 




22 000 




5 x 10 18 


5 x 10 18 


1.0 


2.0 


+25 


+50 


2 


2 


18 000 




3 x 10 18 




1.0 


0.0 







1.5 




18 000 




3 x 10 18 


3 x 10 18 


1.0 


6.5 


+25 


+50 


1.5 


1.5 



our models and our taper criterion, a disk extending in to 
the star's surface along the magnetic plane can protrude no 
further than 3.2i?* from the central plane. For a second case 
for which the disk has an inner radius of 3i?» (i.e., an evac- 
uation of the inner 50% of the disk) , our models produce 
a similarly tapered phase if the disk extends >3.0R» out 
of the magnetic plane. For the cylindrical model we have 
adopted there is a degeneracy between these dimensions - 
we have little leverage in determining where an inner disk 
edge lies. More importantly than this, the implication from 
the cylindrical model representation is that the disk extends 
to surprisingly large distances from the central plane. There 
is no indication that disks of other Bp stars extend out of 
the plane to similar distances. It seems more likely that our 
assumed disk geometry is in conflict with fitting difficulties 
that are aggravated for star-disk systems observed from a 
low inclination. 



Although the concern just raised over our disk geometry 
arises from the combination of high [3 and low i for the C Cas 
case, the resonance line absorption curves of stars having 
high f3 and i values pose their own interpretational dilem- 
mas. The most significant of these is a fact noted by SG01 
for those systems for which the observer's observer's line 
of sight runs through the magnetic plane. In these cases, 
the complete occultation of the star disk in optically thick 
resonance lines presents the physical problem of how supe- 
rions responsible for the line scattering can be distributed 
throughout the cylinder. How can such ions be excited in 
the middle of an optically thick medium? This puzzle can 
be solved if, as predicted in the MHD models of ud-Doula 
and coworkers, the disk is pinched off at the Alfven radius in 
the equator by the confluence of collimated wind streams. 
Exterior to this point the field lines curve out and run par- 
allel to the magnetic equator. Then, an external observer at 
infinity viewing the system in the magnetic plane sees the 



emission from excited superions along this (thin) geometric 
boundary. 

A third difficulty offered for the cylindrical geometry, 
again evident for the C Cas case, is that the absorptions 
should scale differently for optically thick and thin lines as 
the disk moves off the limb of the star. The excess absorp- 
tions of these two groups of lines should scale differently - 
as the area and the volume of the occulting disk, respec- 
tively. In contrast, we find that when the curves in Fig.[T] 
are overplottcd, their morphology is indistinguishable for 
the two types of lines. This difficulty is inherent for any 
"solid" disk geometry. It can be avoided if the disk opac- 
ity is confined to a thin outer region having similar depth 
ranges of formation for these lines. This circumstance is 
permitted if the volumes of excited ions are confined to the 
outer regions of the disk. 

In Figure[4] we suggest a revised geometry in which the 
outer disk edge, representing the wind-disk shock, tapers 
toward the equator and thus takes the general shape of 
an oblate spheroidQ We envisage that the excited "super- 
ions" are formed in a thin outer shell of this structure. The 
column and volume densities can then be made consistent 
if the absorption column lengths, as in Table 1, refer to a 
common shell thickness. For densities predicted by uO02 
and column lengths in Table 1, the shell thickness to t = 
1 is then only of the order 10 3 km. Since the resonance 
lines are formed in a skin around a cool transparent disk, 
there need be no relation between their column densities 
and those inferred from the Al III and CII lines. 

The regions where the resonance line emissions arise 
likewise should be reconsidered, as they cannot be under- 



3 Rather than an oblate spheroid, the circumstellar disk region 
could take on just as well a negative curvature, coming to a cusp 
at the equator. The important point in our picture is not the 
exact shape but that it be closed at the equator and have about 
the height and radii we determine. 
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Fig. 4. The geometry at the times of maximum (</> = 0.5) and minimum (<f> = 0) absorptions of the resonance lines of 
£Cas in panels a) and b), respectively. The observer is to the right. The angles i = 18° and (3 = 79° are taken from 
N03. Low-excitation line absorptions are formed in the optically-thin disk interior. Closed at the outer end, the disk is 
centered around the magnetic equatorial plane and taken to be an oblate spheroid. The outer edge of the disk is heated by 
impacts from the high velocity wind and includes superions responsible for the resonance line absorptions. Net redshifted 
emissions are produced in an "auroral cap" shown as a dotted region over the magnetic poles (at ~11 and ~5 o'clock 
positions). 



stood as the result of the interaction of the wind and the 
edge of a dense disk. Consider first that for the Bp stars 
with dense, extensive disks investigated by SG01, the emis- 
sions can be present at all phases. The SG01 picture was 
that the resonance lines exhibit emission preferentially at 
those phases when the disk is viewed nearly face-on in the 
plane of the sky. A problem with this geometry is the pos- 
itive radial velocity of the emissions. In contrast to SG01, 
we now believe the redshifts do not arise from the impact 
of the wind particles against the disk interface. Such inter- 
actions would produce emissions, but they are centered at 
the rest frame of the star. 

A second problem is that the ratios of volumes com- 
puted from these lines' emissions and excess absorptions 
can be quite diverse among Bp stars. However, they should 
be approximately equal if the contributions in the same 
disk volumes arise from scatterings. These scatterings are 
conserved when integrated over all lines of sight. 

These issues can be addressed by positing that the ab- 
sorptions and most of the emissions of UV resonance lines 
are formed in different places. Our revised geometrical pic- 
ture harkens back to the results of two-dimensional MHD 
models of the emission complex surrounding 6 1 Ori C dis- 
cussed by uO02 and G05. Recall that these models predict 
that infall condensations fall back into the outflowing wind. 
Indeed, the collision of these two structures will produce 
redshifted emission in the external reference frame because 
the center of mass is dominated by the falling blobs. These 
emissions are visible mainly over a narrow range of phases 
because the projected shock areas modulate with viewing 
angle. For 8 1 Ori C the projected radial velocity in the UV 
and X-ray resonance line is observed to be >100 kms -1 
(Smith & Fullerton 2005, G05). 

We now ask what would be outcome of this same pro- 
cess if the wind were weaker? We speculate that the answer 
begins with the fact that because the lower efflux would 
provide less resistance to the infall, the shock front where 
the pressures equalize will lie closer to the star. This shock 
point will then occur where the dense fallback blobs have 



accelerated to a higher velocity, while the wind outflow, in- 
tercepting them from the opposite direction, has only begun 
to accelerate. For a sufficiently weak wind (and similar 77*) 
the shock front should then describe an arc-like "auroral" 
cap not far above the star's magnetic pole with a center- 
of-mass redshift because of the fall back's higher velocity 
and density. We expect that there are many such return- 
ing blobs at any one time, so observations repeated at the 
same phase will see the same level of emission. The pro- 
jected area of the shock viewed by the observer, and hence 
the recorded emission of the resonance lines, would become 
maximum at <j> = 0. (The reverse would be the case at <j) 
= 0.5 when the project area of the cap would be minimal.) 
Clearly, this picture begs to be confirmed by MHD sim- 
ulations of the type conducted for other Bp stars having 
stronger fields. 

The putative auroral shock regions are represented as 
shaded regions at the "11 o'clock" and "5 o'clock" posi- 
tions relative to the star in Fig(4] In this new configuration 
at cf> — these shocks have smaller areas but are denser than 
the the wind-shock interfaces discussed by Babel and SG01. 
Note that because the emissions and absorptions are formed 
mainly in two different regions of the magnetosphere, their 
volumes and column densities need not be equal. In ad- 
dition, one expects that the UV line emission is formed 
in a separate region of the magnetosphere from the X-ray 
emissions. The energies released at the auroral caps can be 
expected to be moderate at most. 

This picture suggests that sometimes one might see 
small amounts of blueshifted emission from the opposite 
hemisphere. As shown by Fig[4] right panel, "4-5 o'clock" 
auroral cap contribution from the aurora over the oppos- 
ing figure would exhibit a much smaller Dopplcr shift and 
projected area than the opposing cap. We can speculate 
that the very high microturbulence taken for conditions in 
a single structure in our previous models, are the result of 
a superposition in velocity of these two shocks (each with 
a smaller £). 
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4. Summary 

Our modeling of the resonance line variations of £ Cas shows 
that the disk occults the star to some degree during the 
entire rotation cycle. This attribute has not been found in 
the quantitative modeling of other magnetic Bp stars whose 
lines we have examined (e.g., SG01, Smith 2003, Smith & 
Fullerton 2005, S06). It is probably more easily found for 
£Cas because of its small inclination angle i. In this ge- 
ometry the disk occultation never realizes the full possible 
range of areal coverages (0-100%). 

As mentioned, the large areal coverage factor found for 
our solutions at both cj) = 0.0 and 0.50 requires that the 
disk either extend nearly to the star's surface, extend to 
large distances from the central magnetic plane, or both. 
In Fig.Hl we noted that most or all the lines of sight to the 
star intersect the disk and thus suffer flux attenuation. 

Consideration of the characteristics of the emission com- 
ponents in ( Cas and the "classical" Bp stars shows that 
they cannot be easily understood with the kinematic de- 
scription adopted by SG01, whether the emission arises 
from a shock at the edge of a cylindical or torus shaped 
disk. For example, for the Bp stars with dense or extensive 
disks investigated by SG01, the emissions can be present at 
all phases, and in all cases with a net redshift. We conjec- 
tured that this might be understood by shocks created by 
the collisions of dense fall back (blobs) with the outflow- 
ing wind. This releases the magnetic disk model from other 
interpretational inconsistencies. These include the frequent 
inequalities between the derived volumes in which the ab- 
sorptions and emission are formed. It also addresses the 
high microturbulcnt velocities that are otherwise necessary 
to fit the line profile variations. 

The analysis of lines of a modestly inclined system like 
£ Cas clarifies these problems and provides cause for a re- 
consideration of the disk and emission region geometries. 
This shows that a final picture of these complicated star- 
disk system requires observations not merely over all view- 
ing angles but of a sample of stars with a variety of rota- 
tional and magnetic inclination angles. 

We thank the referee, Dr. Marc Gagne, for a number of 
comments that have added to the quality of this paper. 
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